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Abstract 

■ As a main example for the superstatistics approach, we study a Brownian particle moving in a 

, d-dimensional inhomogeneous environment with macroscopic temperature fluctuations. We discuss 

I ' the average occupation time of the particle in spatial cells with a given temperature. The Fokker- 

^ ■ Planck equation for this problem becomes a stochastic partial differential equation. We illustrate 

o , 

, our results using experimentally measured time series from hydrodynamic turbulence. 
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I. INTRODUCTION 



Nonextensive statistical mechanics 



mm 



originally developed as an equilibrium 



formalism, has indeed many interesting applications for driven nonequilibrium systems of 
sufficient complexity. Consider e.g. a Brownian particle moving through a changing envi- 
ronment. We may assume that the environment is inhomogeneous and a suitable parameter 
describing the state of this environment, e.g. the inverse temperature P, fluctuates on a 
relatively large spatio-temporal scale. This means that there is a fast dynamics given by 
the velocity of the Brownian particle and a slow one given by temperature changes of the 



environment, 
superstatistics 



le two effectsproduce a superposition of two statistics, or in a short, a 
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The 



stationary probability distributions for an ordinary Brownian particle in a constant envi- 
ronment are certainly Gaussian distributions, but for a fluctuating environment one obtains 
non-Gaussian behaviour with fat tails. These tails can decay e^. with a power law, or as 
a stretched exponential, or in an even more complicated way Which type of tails are 
produced depends on the probability distribution /(/?) of the parameter j3. For the above 
simple example of a generalized Brownian particle, /3 is the fluctuating inverse temperature 
of the environment, but in general (3 can also be an effective friction constant, a changing 
mass parameter, a changing amplitude of Gaussian white noise, the fluctuating energy dis- 
sipation in turbulent flows, or simply a local variance parameter extracted from a signal. 
Recent applications of the superstatistics concept include a variety of physical systerns, such 



as Lagrangian 
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mospheric turbulence 
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and Eulerian turbulence 
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defect turbulence 



33| |. cosmic ray statistics j3^, solar flares 3^, solar wind statistics 



3l|, at- 
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networks j37l . I38| , random matrix theory |39j , and mathematical finance 
If f3 is distributed according to a particular probability distribution, the x^-distribution, 
then the corresponding marginal stationary velocity distributions of the Brownian particle 
obtained by integrating over all (3 aregiven by the generalized canonical distributions of 



aregiven 



nonextensive statistical mechanics For other distributions of the intensive pa- 

rameter (3, one ends up with more general statistics, which contain Tsallis statistics as a 
special case. 

In this paper we first briefiy review the superstatistics concept. We then look at typical 
occupation times of the superstatistical Brownian particle in cells of constant temperature. 
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distinguishing between temporal and spatial temperature fluctuations. In section 4 we em- 
phasize a certain analogy between the step from statistics to superstatistics and the step 
from 1st quantization to 2nd quantization: In both cases the rate equations for probability 

densities are upgraded from deterministic partial differential equations to stochastic ones. 
Finally, in the last section we compare with experimental data in turbulent flows, where 
superstatistical models are a very useful tool to effectively describe the dynamics. 

II. VARIOUS TYPES OF SUPERSTATISTICS 

It is well known that for equilibrium systems of ordinary statistical mechanics the prob- 
ability to observe a state with energy E is given by 

P{E) = ^^P{E)e-P''. (1) 

This formula just describes the ordinary canonical ensemble. e~^^ is the Boltzmann factor, 
p{E) is the density of states and Z{(3) is the normalization constant of p{E)e~^^. For 
superstatistical systems, one generalizes this approach by assuming that /3 is a random 
variable as well. Indeed, a driven nonequilibrium system is often inhomogcncous and consist 
of many spatial cells with different values of (3 in each cell. The cell size is effectively 
determined by the correlation length of the continuously varying /?-field. If we assume that 
each cell reaches local equilibrium very fast, i.e. the associated relaxation time is short, 
then in the long-term run the stationary probability distributions arise as a superposition of 
Boltzmann factors e~^^ weighted with the probabihty density f{/3) to observe some value 
/5 in a randomly chosen cell: 

POO 1 

P{E) = fm^^p{E)e-^^dp (2) 

The simplest example is a Brownian particle of mass m moving through a changing 
environment in d dimensions. For its velocity v one has the local Langevin equation 

il^-^v + aL(t) (3) 

{L{t): (i-dimensional Gaussian white noise, 7: friction constant, a: strength of noise) which 
becomes superstatistical because for a fluctuating environment the parameter '■— 
becomes a random variable as well: it varies from cell to cell on the large spatio-temporal 
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scale T. Of course, for this example E = ^rriip. During the time scale T the local stationary 
distribution in each cell is Gaussian, 

Pirn = [-^Y' e-'^^-^ . (4) 
But the marginal distribution describing the long-time behavior of the particle for t >> T, 

POO 

Piv) = / fiP)pimdP (5) 



exhibits fat tails. The large- |f| tails of the distribution (0) depend on the behaviour of f{(3) 
for (3^0 ^. Many different superstatistical models corresponding to different are 
possible: The function / is determined by the environmental dynamics of the nonequilibrium 
system under consideration. 

Consider, for example, a simple case where there are n independent Gaussian random 
variables Xi, . . . , X„ underlying the dynamics of in an additive way. f3 needs to be positive 
and a positive f3 is obtained by squaring these Gaussian random variables. The probability 
distribution of a random variable that is the sum of squared Gaussian random variables 
(3 = J2i=i -^i is well known in statistics: It is the x^-distribution of degree n, i.e. the 
probability density is given by 

m = ^U^) (6) 



r(f) \2(3o 
where (3o is the average of /3. 

One can now simply do the integration in eq. (0) using eq. (jH)) 3 Q • The result is the 
generalized canonical distribution of nonextensive statistical mechanics, i.e. a g-exponential 
of the form 

p{v) ~ z ^ — (7) 

{1 + p{q - l)lmiP)^ 

with 

q=l + (8) 

n + a 

and 

^ 2-{q-l)d- 

Thus these types of generalized Boltzmann factors occur as stationary states for nonequi- 
librium systems with suitable fluctuations of some intensive parameter (3. 
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Instead of P being a sum of many contributions, for other systems the random variable 
(3 may be generated by multiphcative random processes. In this case one typically ends up 
with a log-normally distributed P [s^, i.e. the probability density is given by 



1 



exp 



-(In^ 



2s2 



(10) 



where /i and s are parameters. Both x^-superstatistics (= Tsallis statistics) and lognormal 
superstatistics are relevant for turbulent systems. The observation is that the former one 
particularly well describes atmospheric turbulence data jS 



33l |. whereas 



;he latter one 
29, m. Note 



gives good results for laboratory turbulence (e.g. Taylor-Couette flow) j2 
that the Reynolds number is controlled for laboratory turbulence, whereas it fluctuates for 
the atmospheric case. 



III. OCCUPATION TIMES IN CELLS OF CONSTANT LOCAL TEMPERATURE 

Let us now discuss some subtleties of the superstatistical modeling approach. An im- 
portant question is what really causes the fluctuations of temperature around the Brownian 
particle. Clearly these are produced by the fact that we consider a driven nonequilibrium 
system in a stationary state which has some external energy input and also some energy 
dissipation (like the sun acting on the earth's atmosphere and producing spatio-temporal 
temperature fluctuations in form of the weather). There are, however, different possibilities: 
Either there could be explicit temporal temperature changes around the Brownian particle, 
or the particle could move through spatial cells of size L'^ which all have different temper- 
atures but no explicit time dependence. The latter case is more a spatial then a temporal 
effect, and corresponds to a frozen random pattern. 

To treat the second case properly, we have to take into account the average occupation 
time of the Brownian particle in a cell of size L'^ where the temperature is approximately 
constant. Locally, in each cell the velocity of our particle is described by the Ornstein- 
Uhlenbeck process and its position by the Ornstein-Uhlenbeck position process [4^. The 
latter process is a Gaussian diffusion process. For each position component Xi(t), i = 1, . . . ,d 
of the Brownian particle one has for large times t 

{x^ = 2Dt, (11) 



5 



where the diffusion constant D is given by the Einstein relation 

Of course, the larger the temperature in a given cell, the faster the particle will diffuse, 
hence it will occupy a given cell of size L*^ for a shorter time on average if the temperature 
is higher in that cell. In fact, the average time t the particle spends in a cell of temperature 
is given by 

= 2Dt ^t=^ = h'm^p. (13) 

Assuming that on average all spatial cells with a given local temperature have the same 
typical size L'^, this means the average occupation time i in a given cell is proportional to 
p. The probability to find the particle in a given cell is proportional to the occupation time, 
hence this means we have to adjust the probability density fspaceiP) describing the spatial 
distribution of /? in the various cells by a factor (3 if we proceed to a temporal description 
ftempW), i-e. if we take into account how much time the particle spends in each cell: 

ftempiP) ~ tfspaceiP) PfspaceiP) (14) 

In general, one has for a given effective energy E the local Boltzmann-Gibbs distributions 

piocim = I3e-^'' (15) 

in each cell (assuming for simplicity a constant density of states). These distributions are 
properly normalized, 

/■oo /"OO 

/ piociE\0)dE^ Pe-f^dE^l. (16) 
Jo Jo 

Taking temporal averages over all possible /3, one obtains the marginal distributions 

/■oo 

p{E) = / (5U,mp{(3)e-P''d(3. (17) 
Jo 

Our previous consideration shows that proceeding to a spatial probability density fspaceiP) 
describing the distribution of /3 in the various spatial cells requires an additional weighting 
by a factor (3, since a typical particle will spend more time in cells with small temperature. 
We thus arrive at 

poo 

p{E) ~ / P'fspUP)e-''''dp. (18) 

^0 
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Our consideration shows that care has to be taken in the interpretation of the probabihty 
density f{(3), namely whether it describes the distribution of jS in spatial cells or in temporal 
time slices. In general, one must also take into account the density of state and the fact 
that for Brownian particles probability densities are usually regarded as a function of the 
velocity v rather than the energy E = ^mip. This corresponds to a simple transformation 
of random variables. 



IV. STOCHASTIC FOKKER-PLANCK EQUATION 

Let us once again emphasize that for a superstatistical Brownian particle one has a 
superposition of two stochastic processes on different time scales: A fast stochastic process 
given by a local Ornstein-Uhlenbeck process and a slow stochastic process given by (3{t). 
The /3-process can a priori be anything. We assume that the parameters a and 7 in eq. (jH} 
are constant for a sufficiently long time scale T, and then change to new values, either by an 
explicit time dependence, or by a change of the environment through which the Brownian 
particle moves. 

For simplicity, let us consider the 1-dimensional case. Locally, in a given cell with constant 
temperature the probability density P{v,t) obeys the Fokker-Planck equation 4^ 

dP divP) I od^P 

with the local stationary solution 



^(^l/5) = ^^exp|-i/3mz;2|. (20) 

In the superstatistical approach, one regards 7(t) and a{t) as stochastic processes as well, 
which evolve on a much larger typical time scale T as compared to v{t). This means that for 
superstatistical systems the above equation is not an ordinary Fokker-Planck equation: It is 
a Fokker-Planck equation whose coefficients 7(t) and a{t) are stochastic processes, with all 
their properties to be specified by a characteristic functional. Hence, in the superstatistical 
approach the equation that determines the evolution of probabilities becomes a stochastic 
partial differential equation — rather than a deterministic one. At the end, one has to perform 
an average over all realizations of 7(t) and a{t). Needless to say that one can generalize all 
this to superstatistical versions of higher-dimensional and nonlinear Fokker-Planck equations 
as well. 
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We can see a certain analogy to the process of second quantization in quantum field the- 
ories. Consider, for example, a Klein-Gordon (or Schrodinger) field. The Klein-Gordon (or 
Schrodinger) equation governs the time evolution of a wave function, and the wave function 
squared essentially gives the quantum mechanical probability density in a first-quantized 
system. Thus the Klein-Gordon (or Schrodinger) equation corresponds to a deterministic 
evolution equation for probability densities. In our case, this role is taken over by the lo- 
cal (deterministic) Fokker-Planck equation. In order to proceed to second quantization in 
quantum mechanics a possible method is to add noise to the classical field equation (e.g. 
the Klein-Gordon equation), obtaining a second-quantized theory via the stochastic quan- 
tization procedure 0,14^. By second quantization, we thus construct a stochastic partial 
differential equation, a Klein-Gordon equation with noise. In a similar way, for superstatisti- 
cal systems we proceed from a deterministic Fokker-Planck equation to a stochastic one, by 
replacing the constant parameters 7 or cr by suitable stochastic processes. Thus proceeding 
to a superstatistical description is analogous to proceeding to a second quantized description 
in quantum mechanics. 

To completely solve the superstatistical system dynamics one has to find all solutions 
of the stochastic Fokker-Planck equation (which are stochastic processes) and then form 
expectations. In full generality, this is quite a heavy task, though some progress has been 
made for very simple dichotomous models Q]. The solution of the stochastic Fokker-Planck 
equation depends on what type of stochastic processes are chosen for ^{t) and cr(t). As shown 
in [2^, one has for example a power-law decay of the marginal velocity correlation function 
if 7 fluctuates and a is constant, and an exponential decay if 7 is constant and a fluctuates, 
choosing the same x^-distribution of /3 = in both cases and assuming T » r, r being 
the local relaxation time in each cell. To construct the best superstatistical model for a given 
physical application, one has to carefully compare with the experimental data, looking at 
both, stationary densities and correlation functions. 



V. COMPARISON WITH TURBULENT VELOCITY FLUCTUATIONS 



A major application of superstatistics is in turbulence modeling 
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I22I, Here essentially the local velocity v of the superstatistical Brownian particle model 
corresponds to a local longitudinal velocity difference u in the turbulent flow, and (3 is a 
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FIG. 1: Probability density ftempl^ as extracted from an experimentally measured time series in 
a turbulent Taylor-Couette flow 30]. Also shown is a lognormal, x^j and inverse distribution 
with the same mean and variance as the experimental data. 



local variance parameter connected with fluctuations in energy dissipation. As said before, 
to construct a suitable superstatistical turbulence model one has to compare carefully with 
various experimentally accessible quantities. For the example of a measured velocity time 
series in a turbulent Taylor-Couette flow, one basically observes the following j^^: The 
probability density ftempiP) is essentially a lognormal distribution (Fig. 1). Correlation 
functions of the measured velocity signal v{t) decay more or less exponentially, but those 
of temporal velocity differences u{t) = v{t + 6) — v{t) on a a given time scale 6 decay 
asymptotically with a power law (Fig. 2). The marginal distribution p{u) exhibits fat tails 
whose flatness decreases with increasing 6. On all scales the measured histograms p{u) are 
well described by lognormal superstatistics as given by eq. © and (jH^), replacing v u. 
One observes a a clear time-scale separation between the process u{t) and the corresponding 



process P{t): f5{t) evolves on a much larger typical time scale, 
increases with increasing Reynolds number, as well as with S 



his time scale separation 



30|. The process P{t) is 



strongly correlated, its correlation function decays with a power law. 
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